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The dynamic world model and its linear perturbations were first studied in Einstein's gravity. In 
the system without pressure the relativistic equations coincide exactly with the later known ones 
in Newton's gravity. Here we prove that, except for the gravitational wave contribution, even 
to the second-order perturbations, equations for the relativistic irrotational zero-pressure fluid in 
a flat Friedmann background coincide exactly with the previously known Newtonian equations. 
Thus, to the second order, we correctly identify the relativistic density and velocity perturbation 
variables, and we expand the range of applicability of the Newtonian medium without pressure to 
all cosmological scales including the super-horizon scale. In the relativistic analyses, however, we do 
not have a relativistic variable which corresponds to the Newtonian potential to the second order. 
Mixed usage of different gauge conditions is useful to make such proofs and to examine the result 
with perspective. We also present the gravitational wave equation to the second order. Since our 
correspondence includes the cosmological constant, our results are relevant to currently favoured 
cosmology. Our result has an important practical implication that one can use the large-scale 
Newtonian numerical simulation more reliably even as the simulation scale approaches near horizon. 
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I. INTRODUCTION 

Despite its algebraic and conceptual complexity in Ein- 
stein's gravity the evolving world model and its linear 
structures were first studied based on Einstein's grav- 
ity in the classic works by Friedmann in 1922 [1] and 
Lifshitz in 1946 [2], respectively. In an interesting se- 
quence, the much simpler and, in hindsight, more intu- 
itive Newtonian studies followed later by Milne in 1934 
[3] and Bonnor in 1957 [4], respectively. In the case with- 
out pressure the Newtonian results coincide exactly with 
the previously derived relativistic ones for both the back- 
ground world model and its first-order (linear) perturba- 
tions. The case with pressure cannot be handled in the 
Newtonian context despite several failed attempts in the 
literature to simulate it especially for the perturbation. 
The situation is still well described by Sachs and Wolfe 
in 1967 [5]: "When these modified equations were per- 
turbed to first order, their solutions did not agree with 
the relativistic results, even qualitatively." In this work, 
we will show an additional continuation of relativistic- 
Newtonian correspondences in the zero-pressure medium 
by proving that the relativistic second-order scalar-type 
perturbations are described by the same equations known 
in Newton's theory. That is, the Newtonian equations 
coincide exactly with the relativistic ones even to the 
second order in perturbations. 

In the relativistic perturbations, due to the covariance 
of field equation we have freedom to fix the spacetime 
coordinate system by choosing some of the metric or 
energy-momentum variables at our disposal: this is often 
called the gauge choice. The original study of Lifshitz 
started by choosing the synchronous gauge which is still 



quite popular in the literature. Other gauge conditions 
were discovered later [6,7]. It is an ironic situation that 
except for the widely used synchronous gauge condition, 
each of other gauge conditions fixes the gauge freedom 
completely. Thus, each has its own unique corresponding 
gauge invariant combination. Notice some common alge- 
braic errors (not in Lifshitz's work though) widespread 
in the literature including many textbooks due to the in- 
complete gauge fixing nature of the synchronous gauge, 
see [8]. 

Although infinitely many gauge conditions are avail- 
able, it has been common in the literature to fix gauge 
conditions from the beginning. The importance of using 
different gauge conditions for different variables and the 
gauge invariance of such variables were shown by Bardeen 
in 1980 [7]. Bardeen's work also showed the importance 
of having access to many different gauge conditions which 
become apparent in his work in 1988 [10]. In this work, 
the importance of having different variables evaluated in 
different gauges (all correspond to unique gauge- invariant 
combinations) will become clear as we extend Bardeen's 
approach to the second-order perturbations. 

Recently, we have presented a second-order perturba- 
tion formulation of the Friedmann world model consid- 
ering quite general situations [11]. We have resolved the 
gauge issue, identifying the variables to use in fixing the 
gauges and constructing gauge- invariant combinations, 
which can be easily extended even to the higher order. 
The basic equations are presented without fixing the tem- 
poral gauge condition thus allowing us to choose or try 
many available gauge conditions later depending on the 
situation: we call this a gauge-ready approach, see Eqs. 
(5)-(ll) below. The Newtonian correspondence to the 
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linear order was made by properly arranging the equa- 
tions using various gauge-invariant variables in [7,12,9]. 
Extending such correspondences to the second order is 
our task in this work. We set c = 1. 



II. BASIC EQUATIONS 

We consider a scalar-type perturbation in the flat 
Friedmann background. We will consider the presence 
of tensor-type perturbation (gravitational waves) in §VI. 
The vector-type perturbation (rotation) is not important 
because it always decays in the expanding phase even to 
the second order, see §VII.E of [11]. Our reason for con- 
sidering the flat background will be explained below Eq. 
(4). As the metric we take 
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ds 2 = -a 2 (1 + 2a) drj 2 - 2a 2 p. a drjdx a 



dx a dx f3 , 



(1) 



which follows from our convention in Eqs. (49), (175), 
and (178) of [11]. Here, a(t) is the scale factor, and 
a, (3, 7 and ip are spacetime dependent perturbed-order 
variables; we take Bardeen's metric convention in [10] 
extended to the second order. A vertical bar indicates 

(3) 

a covariant derivative based on g y a ^ which becomes S a p 
if we take Cartesian coordinates in the flat Friedmann 
background. By taking 7 = 0, which we call the spatial 
C-gauge, the spatial gauge mode is removed completely, 
thus all the remaining variables we are using are spatially 
gauge-invariant to the second order, see §VI.B.2 of [11]. 
In the following we will take 7 = as the spatial gauge 
condition and use \ = a P + fl2 7 which becomes \ = a fi- 
As the energy-momentum tensor we take 



f° = - f x-S l i+-fi X ' a v, a , 
J* 



f$ = 5p6? + 



a) v_ a , 
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-MX' v,p, 



which follows from our convention in Eqs. (84), (175), 
and (178) of [11]; tildes indicate the covariant quantities. 
Here, \x is the background energy density, and S/i, Sp, H 
and v are the perturbed order energy-density, isotropic 
pressure, anisotropic pressure, and the flux, respectively, 
all based on the normal frame vector h a with h a = 0. Al- 
though we are considering a zero-pressure system (thus, 
p = and Sp = = n to the linear order), it is essential 
to keep the perturbed pressure terms Sp and n because 
these do not necessarily vanish to the second order in per- 
turbation depending on the coordinate (gauge) condition 
we choose. This is because in [11] we have evaluated the 
fluid quantities based on the normal- frame n ; we will 
elaborate this point in §111. 

To the background order we have the Friedmann equa- 
tion [1,3,13] 



const. A 
+ 3' 



(3) 



with the energy (mass) density n (g) oc a~ 3 ; A is the 
cosmological constant. To the linear-order perturbations 
we have a second-order differential equation originally 
derived by Lifshitz [2,4] 



S + 2H5- A-kG^S = 0. 



(4) 



An overdot indicates a time derivative based on t (dt = 
adrj) and H = |. The variable a(t) is the scale factor, 

and 5 = ^£ = i£ with /j, (g) and Sfj, (Sg) the background 
and perturbed parts, respectively, of the energy (mass) 
density field. The "const." part is interpreted as the spa- 
tial curvature in Einstein's gravity and the total energy 
in the Newton's gravity [13]. Equation (4) is valid even 
in the presence of the cosmological constant A as well 
as the background curvature. We will include A term 
in the following. In the relativistic context Eq. (4) can 
be derived in the comoving gauge condition; the original 
derivation by Lifshitz is based on the synchronous gauge, 
and in the zero-pressure medium to the linear order the 
synchronous gauge coincides with the comoving gauge: 
further discussion about this point will be made in §IV. 
Although Eq. (4) is also valid with general spatial curva- 
ture the relativistic-Newtonian correspondence is some- 
what ambiguous in the case with curvature, for details 
see §3 of [9]. Thus, we consider the flat background only. 

The perturbed parts of equations to the second order 
are presented in Eqs. (195)-(201) of [11]. In a flat back- 
ground with vanishing background pressure we have 

A 

k - 3Ha + 3(p + -~; 
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(2) k + 2HK-4irG(6ii + 35p)+[3H+-z)a = N 3 , (8) 
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(7) 


A^3, 


(8) 




(9) 


N 5 , 


(10) 
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where the pure quadratic-order terms, Ni, can be read 
from Eqs. (99)-(105) in [11]. A is a Laplacian operator. 
Equation (5) is a definition of k, Eqs. (6)-(9) follow from 
Go, G°, G" - Gfj and Gf 3 - components of Ein- 

stein's equation, respectively, and Eqs. (10), (11) follow 
from f fc ;6 = and f*. b = 0, respectively. To the lin- 
ear order these set of equations without fixing the tem- 
poral gauge was presented by Bardeen in [10]. All our 
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equations include the cosmological constant in the back- 
ground. These equations are presented without fixing 
the temporal gauge condition and using only the spa- 
tially gauge-invariant variables even to the second order; 
our choice of the spatial C-gauge (7 = 0) guarantees such 
invariances of the remaining variables, see §VI.B of [11]. 
As the proper temporal gauge condition we can choose 
any of the following: a = (the synchronous gauge), 
X = (the zero-shear gauge), (5 = (the uniform-density 
gauge), k = (the uniform-expansion gauge), v = (the 
comoving gauge), (p = (the uniform-curvature gauge), 
etc. Except for the synchronous gauge, each of the other 
temporal gauge conditions completely removes the tem- 
poral gauge mode. We can also take linear combinations 
of the above conditions, and choose different gauge con- 
ditions to different order, see §VI.C.2 of [11]. Thus, we 
have infinite number of different temporal gauge choices 
available to each order in perturbations. 

From Eqs. (5)-(ll) we can derive the following set of 
equations expressed using gauge-invariant variables 

a v = ~v x ,a v x , a - ^ (s Pv + |^n^ , (12) 

1 H 

S v - Kv — - (S v v ' a ) - 3—Sp v , (13) 

a y * ' ,a ^ 

A - \ 
k v + 2Hk v - AnGfiSy = (v x a v Xta ) + 12nGSp v , 



A 1 

Kv ~ ~a~ Vx = a ( VxA(fx ~ 2(p * Av x + <Px av x> a ) 



+ -H(2v x Av x +v x a v x , a ) 



1 „ > 3 

a a 
a x + l Px = l P 2 x- A_1 <>x A ^x) 
+3A- 2 V Q V^ (<p x <p x , a p) - 8nGU x , 

AnGfxSy + ^ip x = \h&v x - '5aHH/S.~ 1 \7 a (S ViCt v x 



-^V Q (S v v x . a ) - -A^V" (v x , a Aip x ) , (15) 



(16) 



+ i ( il fx A( Px + ifx^^x^ ) . 



a* 

13- 1 

v x + Hv x - -a x = -^aHv 2 x + 3Hip x v x - —<p x 



(17) 



-lA-^^J + ^^ + f^), (18) 
^x - Ha x + 4nGfiav x = <p x (^p x - ^H(p x ^j , (19) 
<j>v = ^A- X \7 a (v x ' p ip v , a i3 + v x . a A<p v ) . (20) 

Equations (12), (13), (14), and (15) follow from Eqs. (11), 
(10), (8), and (7), respectively, evaluated in the comoving 
gauge. In Eq. (15) we used \v + av x = Xv 1 ^ + a4 ? ' 
and x ( v 9) \v = 0; see Sec. VI.C.2 of [11]. Equation (16) 
follows from Eq. (9) evaluated in the zero-shear gauge. 



Equation (17) follows from Eqs. (6), (7), and using d/i v = 

<5/U — fiav + S[ii q \ ip x = (p — H\ + <Px^ an( i Px^lx = 0- 
Equation (18) follows from Eq. (11) evaluated in the zero- 
shear gauge. Equation (19) follows from Eqs. (5), (7), 
removing k term and evaluating in the zero-shear gauge. 
Equation (20) follows from Eqs. (5), (7), removing k term 
and evaluating in the comoving gauge. In this set of 
equations we located the pure quadratic terms and the 
possible second-order pressure terms on the RHSs. 

Our notation with a perturbed-order variable as a 
subindcx, for example, S v indicates a unique gauge- 
invariant combination of S and v which becomes S un- 
der the comoving gauge condition v = 0. Thus, 5 in 
the comoving gauge is equivalent to a unique gauge- 
invariant combination 5 V . To the linear order we have 
S v = S — a(fi/fx)v. An explicit form of S v to the sec- 
ond order and other gauge-invariant combinations can 
be found in Eqs. (280)- (284) of [11]. As we can construct 
many (in fact, infinitely many) gauge invariant combina- 
tions for S, our notation apparently has the advantage of 
showing explicitly which gauge-invariant combination we 
are considering [14]. 

Here, we briefly discuss a conserved variable to the 
second order. From Eqs. (20), (18), and (16) we have 

^ (a 3 ^)' = - J-^A-iV"^ (<p v , a <p v ,(,) • (21) 



(14) To the linear order we have 

Ifiy = C(X). 



(22) 



Thus, ip v remains constant in time. In the large-scale 
limit (super- horizon scale), ignoring the quadratic-order 
spatial gradient terms, Eq. (22) remains valid even to 
the second order; for more general proof considering the 
pressure term see [15,11]. 



III. ISSUE OF PRESSURE 

Now, we discuss the role of pressure terms in a medium 
without pressure. From Eqs. (233), (235) of [11] we notice 
that the gauge (coordinate) transformation to the second 
order causes pressure (both isotropic and anisotropic) 
terms to appear even in the case without pressure orig- 
inally (physically). Such a complication occurs because 
our fluid quantities introduced in [11] are based on the 
normal-frame four-vector n a which differs from the fluid 
four- vector u a . In [11] we have presented the fluid quan- 
tities based on u a separately as well, see Eqs. (87), (88) of 
[11]; by using these equations we can translate fluid quan- 
tities in the normal frame to the ones in the fluid frame, 
and vice versa; the gauge transformation properties of 
the fluid quantities in the fluid frame are presented in 
Eq. (238) of [11]. The isotropic and anisotropic pressures 
are gauge (coordinate) dependent quantities. To the lin- 
ear order in the Friedmann background the anisotropic 
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pressure is gauge invariant and the perturbed isotropic 
pressure depends on the coordinate only if we have non- 
vanishing (and time varying) background pressure. In 
the normal-frame, the pure coordinate transformation to 
the second and higher orders will cause both pressures 
(i.e., isotropic and anisotropic pressure like terms in the 
energy-momentum tensor) generated even in the case of 
vanishing pressures to the background and to the lin- 
ear order, see Eq. (233) of [11]; the frame dependence of 
fluid quantities was studied in [16]. This complication 
does not occur for the fluid quantities based on the fluid 
frame-vector u a ; see Eq. (238) in [11]. 

For vanishing pressure terms in the background and 
first-order perturbations we have the following gauge- 
invariant combinations of pressure terms (based on n a ) 
[17] 



four- vector h a . In the normal-frame information about 
the fluid motion is present in the flux four-vector q a with 
q a h a = 0. In the energy frame, which takes vanishing 
flux q a = as the frame condition, the comoving gauge 
condition takes u a = for the fluid four-vector; here, 
we ignore the vector-type perturbation. Since u a = it 
coincides with the normal frame vector. Now, in the nor- 
mal frame, which takes n a = as the frame condition, 
the comoving gauge condition without rotation implies 
q a = 0. Thus, as long as we take the comoving gauge 
without rotation, in either frame we have q a = and 
u a = = h a ; i.e., the fluid four- vector coincides with 
the normal four- vector. 



IV. A PROOF 



6p v = 5p- -p.v a v t<x , 
n„ = n - |/«i 2 A- 2 V a V /S f v, a v. 



(23) 



From this we notice that the gauge- invariant combination 
Sp v is the same as 5p in the comoving gauge. Evaluating 
Eq. (23) in the zero-shear gauge (x = 0) and using v x = 
v — -v to the linear order, we have 



Sp x = Spv + -fJ,v x ' a v XjCt , 



u x = u v 



(24) 

As the definition of fluid without pressure we set the pres- 
sure terms in the comoving gauge equal to be zero, thus 



Sp v = = IT,,, 



(25) 



which are gauge-invariant (and physical) zero-pressure 
conditions. Thus, 

Sp x = ^ v x av x-,c<, 

II X = ^VA^V-V 3 {v x<a v x ,p - ^i> x 'X,7) • 

(26) 

We set the pressure terms using Eqs. (25), (26). Thus, 
for fluid quantities based on the normal-frame, in the 
gauge other than the comoving gauge the physical zero- 
pressure condition implies presence of pressure terms in 
the definition of the energy-momentum tensor. 

In the comoving gauge without rotation the two 
frames, u a and n a , coincide. The normal frame h a has 
n a = 0. The fluid quantities are ordinarily defined in the 
fluid (u a ) frame which differs in general from the normal 



Now, we come to our main point proving the rela- 
tivistic -Newtonian correspondence to the second order. 
Combining Eqs. (13), (14) we can derive [18] 



5 V + 2H5 V - 4ttGhS v 

1 d r fx *\ 



2^ IV" 



(27) 



Equations (13), (14), (17), and (27) can be compared 
with the Newtonian perturbation equations. 

The mass conservation, the momentum conservation, 
and the Poisson's equation in Newtonian context give [19] 



• 1 

5 + -V-U 

a 

u + Hu + -V(5$ 

a 

From these we have 
<5 + 2HS - AttGqS 

— ps^-w: 



-iv-(«u), 

a 
1 

— u- Vu, 

a 

A-kGqS. 



(28) 
(29) 
(30) 



r V • (u • Vu) . 



(31) 



In the Newtonian context Eqs. (28)-(31) are valid to 
fully nonlinear order; i.e., the zero-pressure Newtonian 
fluid equations are exact in quadratic order nonlincar- 
ity Equation (31) has been analysed extensively in the 
Newtonian context, see [20,21]. 

To the linear order it is well known that 5 V , — Vw x and 
— tp x (or a x ) correspond to a density perturbation (5 = 
^ with q = g + Sg and g the mass density), a velocity 
perturbation (u) and a perturbation of the gravitational 
potential respectively. [7,12,9]. To the linear order 

we may identify [9] 



S = 6 V , (5$ 
u = — Vu v , 



-fx 



a 



X' 



1 A 

V • U = — V v = K v . 



(32) 
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As we identify 

6 V = S, k v = — -V • u, (33) 
a 

to the second order, Eq. (27) coincides exactly with Eq. 
(31). Equation (13) becomes 

*„ + -V-u=--V-(<y„u), (34) 
a a 

which coincides with Eq. (28). Equation (14) gives 
V • (u + Hu) + A-nGiJLa5 v = ~V-(u-Vu), (35) 

which also follows from Eqs. (29), (30) in the Newtonian 
context. This completes our proof of the correspondence. 
Such identifications of density and velocity perturbations 
imply that we cannot identify — ip x (or a x ) with 5$ to 
the second order. This conclusion follows from a close 
examination of Eqs. (12)-(20). In fact, using the intrinsic 
three-space curvature in Eq. (265) of [11] 

R {h) = ^| [-2A</? + 8^A<^ + 3<^V, a ] , (36) 
Eq. (17) becomes 

AttGh5 v - ^R[ h) = \hAv 2 x - 3aHHA~ 1 V a (5 v , a v x ) , 

(37) 

which still differs from the Newtonian Poisson's equation. 
Thus, we conclude that we do not have a relativistic vari- 
able which corresponds to the Newtonian potential to the 
second order. Apparently, it is essentially important to 
employ mixed gauge conditions, i.e., take different gauge 
conditions for different variables, to make correspondence 
with the Newtonian system: in this way, correct identi- 
fications of (gauge-invariant) variables are important to 
show the relativistic-Newtonian correspondence. 

At this point, let us clarify the meaning of the quanti- 
ties involved in Eqs. (32), (33). Variables a, \ an d <P are 
defined in the metric in Eq. (1). Variables x an d <P can be 
further identified as the perturbed shear and perturbed 
three-space curvature of the normal hypersurface, respec- 
tively. From Eq. (36) we find that the intrinsic scalar 
curvature i?^) vanishes for tp = 0. From Eq. (264) of 
[11] we find that the tracefree part of the extrinsic curva- 
ture tensor K a p (equivalently, shear tensor of the normal 
frame vector with a minus sign) vanishes for \ = 0, The 
variable k can be interpreted as the perturbed expansion 
with a minus sign. From Eqs. (57), (99), and (179) of 
[11] we have K = —3H + n where K is a trace of the 
extrinsic curvature tensor K a p (equivalently, the expan- 
sion scalar, 6 = fi a . al with a minus sign). Variables 8 and 
v are defined in Eq. (2) and can be interpreted as the 
perturbed energy-density (5 = ^ with ft = /i + 5[i) and 
the flux of the normal-frame, respectively. In the normal 



frame, from Eqs. (4), (76), and (175) of [11] we have the 
flux vector becomes J a = — n^T^ = —a^iv_ a . 

Here we discuss the relation between the comoving and 
the synchronous gauge to the second order. Equation 
(12) shows that a v , which is the same as a in the comov- 
ing gauge (v = 0), docs not vanish to the second order. 
This means that the comoving gauge does not imply our 
synchronous gauge to the second order in a zero-pressure 
medium. At this point it is important to remember that 
we already have fixed the spatial gauge condition using 
7 = 0. The original synchronous gauge used by Lifshitz 
fixes Sgoo = = 8go a , thus a = for the temporal gauge 
and (3 = for the spatial gauge condition. We prefer to 
fix 7 = (spatial C-gauge) as the spatial gauge condi- 
tion instead of (3 = (spatial -B-gauge) because the latter 
condition fails to fix the spatial gauge degree of freedom 
completely whereas the first one fixes it completely; this 
is true even to the second order, and in fact to all or- 
ders, in perturbations, see §VI.B.2 and VI. C of [11]. We 
can show that the comoving temporal gauge (v = 0) to- 
gether with spatial i?-gauge (/? = 0) implies a = even 
to the second order, for a proof see [22] . By imposing the 
comoving (v = 0) and the synchronous (a = 0) gauge 
conditions simultaneously, Kasai [23] has derived a dif- 
ferent equation compared with ours: such a redundant 
choice is allowed as one takes /3 — as the spatial gauge 
condition. However, in that gauge condition (the spatial 
B-gauge) the spatial gauge-mode is incompletely fixed, 
and the comparison with the Newtonian analyses is not 
available. 



V. FULLY NONLINEAR EQUATIONS 

By extending our comoving gauge condition to be valid 
to all orders we can formally derive the completely non- 
linear equations for the density and velocity perturba- 
tions. We will present two methods to reach such non- 
linear equations. These are based on the ADM (3 + 1) 
equations and the covariant (1+3) equations summarised 
in §11. A and II. B, respectively, of [11]. With the hindsight 
gained from our second-order perturbations in previous 
sections, it is best to take the comoving gauge condi- 
tion to all orders. In the normal- frame context, only 
the comoving gauge allows the zero-pressure conditions 
to be, by definition, vanishing pressure terms to all or- 
ders. To the second order, all the equations we need to 
derive Eqs. (27), (34), and (35) are Eqs. (12), (13), and 
(14) which follow from Eqs. (8), (10), and (11); these are 
the Raychaudhury, the energy conservation and the mo- 
mentum conservation equations, respectively. We have 
presented a redundant set of equations in (12)-(20) in or- 
der to show the relativistic-Newtonian correspondences 
with some perspective. 

The complete set of ADM (3 + 1) equations is pre- 
sented in Eqs. (8)-(13) of [11], see [24] for original studies. 
We only need Eqs. (10), (12), and (13) of [11] which are 
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the trace of ADM propagation equation, and the energy 
and momentum conservation equations, respectively. We 
take the comoving gauge condition to all orders which 
makes the flux four-vector to vanish, i.e., J a = 0; here 
we assume vanishing vector-type perturbation, thus ir- 
rotational, which could contribute to J a . Under such 
conditions the zero-pressure conditions (in our normal 
frame) imply S = = S a p to all orders; S and S a p are 
the trace and tracefree parts, respectively, of the spatial 
part of energy-momentum tensor. Equation (13) of [11] 
gives 



0, 



(38) 



where N is defined as g oa = —N 2 . Thus, we may set 
N = a(t) to all orders. In this case we have, for example, 
E = E^N- 1 . Now, Eqs. (12), (10) of [11] become 



E - KE = 0, 



(39) 



K - -K 2 - K aP K aP - AnGE + A = 0, (40) 

o 

where E = E — E_ a N a N~ 1 , etc.; E is the energy density 
based on normal frame vector, and K and K a p are the 
trace and tracefree parts, respectively, of the extrinsic 
curvature; N a is defined as goa = N a . The spatial indices 
in ADM formulation are based on the spatial metric h a p 
defined as h a p = g a p. By combining these equations we 
have 



K af3 K afj - AttGE + A = 0. (41) 



Notice again that Eqs. (39)-(41) are valid to all orders, 
i.e., these equations are fully nonlinear. From Eqs. (39)- 
(41), using 




E = Sfj,, 



(42) 



and the quantities presented in [11] we can easily derive 
Eqs. (34), (35), and (27), respectively; see the next sec- 
tion. 

The complete set of covariant (1 + 3) equations is pre- 
sented in Eqs. (26)-(37) of [11]; see [25] for original stud- 
ies. We only need Eqs. (26)-(28) of [11] which are the 
energy and momentum conservations and the Raychaud- 
hury equation, respectively. We take the energy-frame 
which sets the energy flux term to vanish, i.e., q a = 0. 
In this frame the frame four-vector u a is the fluid four- 
vector. The zero-pressure conditions imply p = = TT ab 
to all orders; ir ab is the covariant anisotropic stress based 
on u a . Equation (27) of [11] gives vanishing acceleration 
vector, i.e., a a = u a:b u b = to all orders. Thus, Eqs. 
(26), (28) of [11] become 

jl + 0= 0, (43) 
9+^9 2 + ~o ah ~o ab - u ab £j ab + inGfl - A = 0, (44) 



where 9 = u a . a is an expansion scalar, and d ab is the shear 
tensor. An overdot with_tilde is a covariant derivative 
along the u a vector, e.g., fi = jl,aU a . By combining these 
equations we have 




~a ah ~a ah + £u ab u ab - AvGfi + A = 0. 

(45) 



Notice that Eqs. (43)-(45) are valid to all orders, i.e., 
these equations are fully nonlinear. More general equa- 
tion in a fully covariant form considering the general pres- 
sure terms can be found in Eq. (88) of [26]. 

We take the comoving gauge condition to all orders 
which makes the space part of four-vector with low in- 
dex to vanish, i.e., u a = 0; here we also assume vanishing 
vector-type perturbation, thus irrotational, which could 
contribute to u a . As our gauge condition (and the irrota- 
tional condition) implies u a = 0, the frame vector is the 
same as the normal frame, thus u a = n a . In such a case 
we have vanishing rotation of the u a flow, thus u) a b = 0. 
From Eqs. (43)- (45), using 



ft = n + Sfi, 



(46) 



and the quantities presented in [11] we can easily derive 
Eqs. (34), (35), and (27), respectively. A derivation based 
on the covariant equations is presented in [27]. 

Afterall, the ADM equations (39)-(41) are the same as 
the covariant equations (43)-(45), expressed in different 
forms. We can derive the ADM equations by rewriting 
the covariant equations in the normal frame vector. Since 
our comoving gauge condition with irrotational condition 
makes u a = 0, the frame vector is the same as the normal 
frame vector h a . By direct calculations, using the quan- 
tities presented in Eqs. (2)-(6), and (14)-(16) of [11], we 
can show that 

/i = £, 9 = -k, a ab a ab = K a ?K a p, 

jl = E, 9 = -K. (47) 

Using this we can show that Eqs. (43)- (45) give Eqs. 
(39)-(41); these equations are valid considering general 
background curvature and the tensor-type perturbation 
(gravitational waves) to all orders. 



VI. ANOTHER DERIVATION INCLUDING THE 
GRAVITATIONAL WAVES 

Since Eqs. (34), (35) are our main results allowing us 
to conclude about the relativistic-Newtonian correspon- 
dence, in the following we will derive these equations in 
some detail again directly from the fully nonlinear equa- 
tions in §V. Now, we include the gravitational wave con- 
tribution. The metric becomes 
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ds 2 = -a 2 (1 + 2a) drj 2 - 2a X , a dr)dx a 



+a' 



(1 + 2<p) S af3 + 2<7<g dx a dxP, (48) 



where C^g is the transverse and tracefree gravitational 

(3) 

waves; its indices are based on g y a p. We work in the 

temporal comoving gauge. Thus, C^g is also evaluated in 
the comoving gauge, and equivalent to a gauge-invariant 
combination C 



(*) 



y af3v' 

We introduce 



E = n + 8n, K 



(49) 



see Eqs. (45), (72), (178), and (179) of [11]. We have 
E = E- E !Ct N a N~ 1 =fi + 8ii+ \8ti, aX ' a , 



K = K — K a N a N~ 1 = -3 



+ k+\n, aX a . (50) 



£ + 3- ] (1 + 8 V ) + 8 V + -V • u = --V • (8 v u) , (55) 

3- + 47rG/i - A + -V ■ ( u + -u ) + 4ttGli8 v 
a a \ a J 

= -^V (u • Vu) - (H Ua>/3 + C«) . (56) 

The perturbed parts give Eqs. (34), (35) with additional 
contributions from the gravitational waves in Eq. (35), 
thus in Eq. (27) as well. 

Therefore, in the presence of the tensor-type perturba- 
tion we have 

5„ + ^V-u = -^V-(5„u), (57) 
l -V ■ ( u + -u ] + A-KGfi8 v = — • (u • Vu) 



(58) 



In setting N = a we already have used the comoving 
gauge condition. Since we take the comoving gauge we 
often ignore the subindex v which indicates the comov- 
ing gauge choice (equivalently the unique corresponding 
gauge- invariant combination between the variable and v); 
for example, our 8 is the same as a gauge-invariant com- 
bination 8 V which is the same as 8 in the comoving gauge 
setting v = 0. Using Eqs. (55), (57), and (175) of [11] we 
can show 



k a0 k af3 = 



i 



MP 



X,a|/3 " ^ (Ax) 2 



Equations (39), (40) become 



a 



3-)(1 + 8)+8-k = k8 



a 



(51) 



(52) 



3- + inGn - A - k - 2-k + 4nGu8 
a a 



1 



1 



a 1 3 a 4 



X^ l/3 X, Q |/3 - i (Ax)' 



(53) 



Now, we have to relate x(= Xu) to our notation. Ap- 
parently, we need x only to the linear order. To the linear 
order the G° a -component of Einstein equation in Eq. (15) 
gives = -k v = £V-u; we have \v = X~av = -av x 
to the linear order. As our u is of the potential type, i.e., 
of the form u = u Q , we have 



u 



-Vx„, 
a 



(54) 



to the linear order. Thus, we have 



thus 



CL 1 d 

8 V + 2-8 v - 4nGfi8 v = - — — [aV ■ (8 v u) 
a a (jTj 



+ J_V . (u • Vu) + C^ aSi (^u a , + C, 



(59) 



The presence of linear-order gravitational waves can gen- 
erate the second-order scalar-type perturbation by gen- 
erating the shear terms. Here, we note the behaviour 
of the gravitational waves in the linear regime. To the 
linear order the gravitational waves are described by the 
well known wave equation [2] 



al3 a af} ~ a 2 al3 ~ 



(60) 



In the super-horizon scale the non-transient mode of C t 



it) 

a/3 



remains constant, thus C a p = 0, and in the sub-horizon 

scale, the oscillatory C^p redshifts away, thus C^p oc a -1 . 
Thus, we anticipate that the contribution of gravitational 
waves to the scalar-type perturbation is not significant to 
the second order. 

To the second order the equation for tensor-type per- 
turbation (gravitational waves) can be derived from Eqs. 
(103), (210) of [11]. Since we are ignoring the vector-type 
perturbation from Eqs. (211), (199) of [11] we have 



K 



-\ (v a V, - ^$a) A~ 2 V^V s N 4lS , (61) 

where we assumed a fiat background and set anisotropic 
stress to be zero. From Eq. (103) of [11] to the second 
order we have 
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N Aa p = ^< 



^2 + <P,(aX,P)) + 



1 x 

+ -lX' 7 | ce X, 7 | / 3 + -J («X,a|/3 ~ 4 W,«|^ ~ 3^, a ^) 



2^ |a C«-2A^$-4pAC$ 



\ ~f(a\fj)d afj\jd ~fd\af3 J 

+ 2 C ^(lv., |i + C S )]}' + i 



'/3|7 



4 X' 7 ' X,7|5 



+ - 



(62) 



In Eq. (62) we have ignored a and ip terms which are 
already quadratic order in the comoving gauge, see Eqs. 
(12), (20). Since we are in the comoving gauge, we have 

X = Xv, V = <Pv, « = k v and C^l = C^ jv . Apparently, 
we need \vi K v and (p v to the linear order. We have 
k v = — -V • u and u = j;Vxt>. For <p v we have 



= <p — aHv = ip x — aHv x , 



(63) 



where we have (p x = — <5<I> and u = — Vw x in Eq. (32). 
Using these identifications we can express the scalar-type 
perturbation variables in Eq. (62) in terms of the New- 
tonian variables. 

Equations (57), (58), and (61) provide a complete set 
describing the scalar- and tensor-type perturbations to 
the second order in the flat Friedmann background. From 
these equations we can see that the linear-order scalar- 
type (tensor-type) perturbation works as a source for the 
tensor-type (scalar-type) perturbation to the second or- 
der. Such effects and the presence of the gravitational 
waves are purely general relativistic ones. 



VII. DISCUSSION 

We have shown that to the second order, ignoring 
the gravitational wave contribution, the zero-pressure 



relativistic cosmological perturbation equations can be 
exactly identified with the known equations in Newto- 
nian system, compare Eqs. (57)-(59) with Eqs. (28)-(31). 
More precisely, the relativistic equations can be identi- 
fied with the continuity equation and the divergence of 
the Euler equation replacing the Newtonian gravitational 
potential using Poisson's equation. In order to achieve 
such a correspondence we need correct identification of 
gauge-invariant density and velocity perturbation vari- 
ables as in Eqs. (32), (33). It is important to notice that 
we have avoided using the potential-like variable in our 
identification. In fact, we showed that we do not have a 
relativistic variable which corresponds to the Newtonian 
potential to the second order. This is understandable be- 
cause the gravitational potential introduced in Poisson's 
equation reveals the action-at-a-distance nature and the 
static nature of Newton's gravity theory compared with 
the relativistic gravity. 

As a consequence, to the second order, the Newto- 
nian hydrodynamic equations (31), (34), and (35) re- 
main valid in all cosmological scales including the super- 
horizon scale. Although showing the equivalence of the 
zero-pressure relativistic scalar-type perturbation to the 
Newtonian ones to the second order, may not be en- 
tirely surprising it should not be so obvious either. It 
might be as well that our relativistic results give rela- 
tivistic correction terms appearing to the second order 
which become important as we approach and go beyond 
the horizon scale. Our results show that there are no 
such correction terms appearing to the second order, and 
the correspondence is exact to that order. A comple- 
mentary result, showing the relativistic-Newtonian corre- 
spondence in the Newtonian limit of the post-Newtonian 
approach, can be found in [28]. In fact, the Newtonian 
hydrodynamic equations appear naturally as the zeroth- 
order post-Newtonian limit [29]. 

We note that although we assumed a flat background, 
our equations are valid with the cosmological constant. 
Thus, these are compatible with current observations of 
the large-scale structure and the cosmic microwave back- 
ground radiation which favour near flat Friedmann world 
model with non- vanishing A [30]. As we consider a flat 
background the ordinary Fourier analysis can be used to 
study the mode-couplings as in the Newtonian case in 
[21]. Our result also may have the following important 
practical cosmological implication. As we have proved 
that the Newtonian hydrodynamic equations are valid in 
all cosmological scales to the second order, our result has 
an important cosmological implication that large-scale 
Newtonian numerical simulation can be used more reli- 
ably in the general relativistic context even as the simula- 
tion scale approaches near (and goes beyond) the horizon 
scale. 

At this point, it is important to be reminded that we 
showed the relativistic-Newtonian correspondence for the 
density and velocity perturbations, but not for the gravi- 
tational potential. Therefore, although our result assures 
that one can trust cold dark matter simulations at all 
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scales for the density and velocity fields, it does not imply 
that one can trust the Newtonian simulations for effects 
involving the gravitational potential, like the weak grav- 
itational lensing effects. Indeed, in order to handle the 
lensing effects properly we often require an extra factor 
of two which comes from the post-Newtonian effects.* 

Since the Newtonian system is exact to the second or- 
der in nonlinearity, besides the gravitational wave con- 
tribution to the second and higher order, any nonvan- 
ishing third and higher order perturbation terms in the 
relativistic analysis can be regarded as the pure relativis- 
tic corrections. Expanding the fully nonlinear equations 
in (43)-(45) or (39)-(41) to third and higher order will 
give the potential correction terms. Our recent investi- 
gation of this important open question shows that to the 
third order there occur pure relativistic correction terms 
which are of ip v -order higher [31]. Thus, the corrections 
are independent of the horizon and are small; see the 
accompanying contribution in [31]. 

In this work we have considered an irrotational sin- 
gle component dust in the flat background. Extending 
any of these assumptions could lead to situations which 
deserve further attention. First, it would be interesting 
to see up to what point the correspondence between the 
two theories can be extended in the non-flat case. In this 
way we can identify possible relativistic effects caused 
by the non-flat nature of the background. Second, in 
this work we have ignored the vector-type perturbation 
because it simply decays in the expanding phase. This 
has to do with considering only the longitudinal part of 
u in Eqs. (35), (58). It would be interesting to include 
the rotational mode to see the similarity and difference 
between the two gravity theories. As the realistic New- 
tonian simulations include the whole u vector as the per- 
turbed velocity it would be practically important to see 
the role of relativistic vector-type perturbation to the 
second order, and to determine whether the relativistic 
effect could be important. Third, the usual cosmolog- 
ical simulations include the cold dark matter together 
with the baryon, thus a system with two components. 
Thus, the relativistic nonlinear perturbations of the zero- 
pressure but multi-component system would be interest- 
ing subject in practice. It is, a priori, unclear whether the 
relativistic-Newtonian correspondence would continue in 
such a multi-component case. In the second and the third 
subjects, the comoving gauge issue should be applied 
with care. Fourth, the presence of substantial amount of 
pressures (both isotropic and anisotropic) would lead to 
relativistic corrections. Even in the linear perturbation 
the pressure terms cause new relativistic correction terms 
which arc not present in the Newtonian system. Thus, 
including the pressure terms in relativistic second-order 
perturbation is interesting because most of the terms will 



*We thank the anonymous referee for suggesting this point. 



be pure relativistic corrections. Such a formulation would 
be practically interesting because we anticipate presence 
of strong pressure in the early universe. All these four 
subjects are left for future studies. 
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